Controlling dipole-dipole frequency shifts in a lattice-based optical atomic clock 
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Motivated by the ideas of using cold alkaline earth atoms trapped in an optical lattice for real- 
ization of optical atomic clocks, we investigate theoretically the perturbative effects of atom-atom 
interactions on a clock transition frequency. These interactions are mediated by the dipole fields as- 
sociated with the optically excited atoms. We predict resonance-like features in the frequency shifts 
when constructive interference among atomic dipoles occur. We theoretically demonstrate that by 
fine-tuning the coherent dipole-dipole couplings in appropriately designed lattice geometries, the 
, undesirable frequency shifts can be greatly suppressed. 
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The development of increasingly accurate atomic clocks has led to many advances in technology and tests of funda- 
mental physics. In the search for the next generation of clocks and frequency standards, there has been considerable 
interest in using alkaline earth species because of their narrow intercombination lines in the optical spectrum 0. In 
\ order to achieve a high level of short-term stability and long-term reproducibility and accuracy on the clock transition, 
it is desirable to have a large number of cold atoms located in a well-characterized trap for an improved signal-to-noise 
qq j ratio {S/N) and reduced systematic errors associated with atomic motion. Single ion-based systems do effectively 
. eliminate Doppler and other motion-related systematic errors when the single ions are confined in the Lamb-Dickc 
regime 0, although the achievable S/N is limited by single quantum absorbers. For neutral atoms it is important 
that changes in the level structure due to the trapping potential do not alter the relevant clock transition frequency. 
Such a scheme has been proposed by trapping alkaline earth atoms in three-dimensional optical lattices tuned to a 
"magic" wavelength where the relevant states for the clock transition experience exactly the same level shift Q . The 
1 S (F = 9/2) - 3 P (F = 9/2) forbidden transition (A = 700 nm) in 87 Sr Q is in particular a promising candidate for 
' a lattice-based optical clock transition because of the long lifetime of the excited state (~160 s) and the insensitivity 
of the J = states to thepolarization state of the trapping light. Already there have been efforts towards the cooling 
and trapping of 87 Sr 0. an d recently this transition was directly observed and measured for the first time 0. 

Calcium, another alkaline earth atom that has been studied extensively as a frequency standard 0, , may be a 

candidate for optical lattice clocks as well. 

In the case of A independent atoms, one benefits from a y/~N improvement in S/N in spectroscopy. However, 
r* ■ atoms trapped in an optical lattice can interact with each other and cannot truly be considered independent. Each 
optically excited atom represents essentially a point dipole whose radiated electromagnetic field can affect other atoms. 
These atom- atom interactions can manifest themselves as shifts in the observed transition frequencies. Because of the 
spatial ordering of atoms in a lattice and the potentially high atomic density, it is possible that such interactions may 
produce very large frequency shifts. One might expect then that dipole-dipole interactions can be much more severe 
here than in, for example, atomic fountains, and thus could place serious limits on the accuracy of an optical lattice 
clock if not accounted for. On the other hand, it might be possible to design lattice geometries where this shift is 
reduced or cancelled. Although the trapping lasers are constrained to operate at the "magic" wavelength, the lattice 
geometry can be altered by changing the relative orientations of the trapping beams, whose degrees of freedom can 
be characterized by a set of variables {a}. 

In this paper, we investigate theoretically the dipole-dipole interaction-induced shifts in the clock transition fre- 
quency recovered by Ramsey spectroscopy. We show that by varying the lattice geometry we can quantitatively 
control the clock frequency shift and even reduce the shift to zero. In particular, we give an analytical equation that 
can be solved giving configurations {ao} where constructive interference causes the line shift to be very large. In 
these "bad" lattice configurations, the magnitude of the line shift scales approximately like A 2 / 3 . Quite generally 
we propose that by tuning the parameter space {a} to lie in between two of these bad configurations, one can find 
"good" configurations where the shift is cancelled. The mechanism of cancellation is associated with the destructive 
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interference of contributions to the shift from different atoms in the lattice. 

It is important to emphasize that the present mechanism and theoretical treatment differ considerably from the 
conventional approaches used to treat dipole-dipole line broadening and shifts. In the case of atoms in a spatially- 
ordered lattice geometry, long-range effects are important, and the usual methods involving binary collisions of nearest 
neighbors ^(| are n °t applicable. These long-range effects include, in particular, interference of the far-field dipole 
radiation produced by the excited atoms, a phenomenon similar to Bragg scattering in a crystal. 

This paper is organized as follows. In Sec. ^ we derive equations describing the evolution of an atomic system 
with dipole-dipole interactions. These equations are derived assuming that the atoms are in the Lamb-Dicke regime, 
with one atom or less per lattice site. In Sec. IIIII we give a brief review of Ramsey spectroscopy and solve for the 
dipole-dipole induced line shift using perturbation theory. We find that the shift can be qualitatively understood 
in terms of the classical interaction energies between oscillating dipolcs. There is a contribution to the shift that 
is zeroth order in the interrogation time t, which is due to imperfections in the Ramsey pulses. Even with perfect 
pulses, however, one finds a shift that is first order in t that results from spontaneous decay of the atoms. Sec. IIVI 
discusses how our result for the line shift can be generalized for systems with imperfect filling of the lattice sites and 
for multilevel atoms. In the case of imperfect filling, one can calculate the mean value of the frequency shift as well as 
some nonzero variance, due to the uncertainty of how the lattice is filled. In Sec.[V] we derive an equation that can be 
solved giving lattice configurations where the shift is large due to constructive interference. We derive an approximate 
scaling law for the shift in these "bad" configurations and discuss how the line shift can be reduced by choosing an 
appropriate lattice design. In Sec. I VII we demonstrate these results numerically for one specific lattice configuration. 



II. EQUATIONS OF MOTION 



To treat the problem of interacting atoms in a lattice, we consider N two-level atoms in the Lamb-Dicke limit with 
polarizability along the z-axis. A simple model of the system consists of treating the atoms as point dipoles, and we 
further assume that there is one or less atom per lattice site. This corresponds to a Mott-insulator state for bosons or 
a normal state for fcrmions. In principle, to solve exactly the problem of interacting atoms one would start from the 
full atom-field Hamiltonian and take into account not only all the atomic degrees of freedom but the continuum of 
electromagnetic field modes. To simplify the theoretical treatment, we effectively eliminate the field in the standard 
way using the Born-Markov approximation (see Appendix). This is valid provided that the atomic system evolves 
slowly on timescales of the correlation time r c , which is of the order L/c where L is the linear size of the system and 
c is the speed of light. As a result of eliminating the field, one finds an effective equation of motion for the density 
matrix p of the atomic system. Atom-atom interactions then appear through an effective Hamiltonian H c s as well as 
through a non-Hermitian operator C: 

^ = l [Ho + H cS ,p}+C[p}. (I) 
ot in 

Here, Hq is the atomic Hamiltonian for a non-interacting system. Writing out all the terms in detail, 
dp l r No , 1 iT \^ n \ r + - 1 

a a^b 

~ £ (KV.p) - 2<r 6 -^+) - \ $>(p - <p<) , (2) 



a, b 

where 

/(v) = | ^sin 2 0- + (3 cos- ti - I ) 

/x 3 / . 9 „cos v ,„ 9 ( cosw sinu\ \ . . 

g(y) = -- {sin 2 6— + (3 cos 2 6- I) + —JJ , (3) 

and 9 is the angle that v makes with the z-axis. 

The first term on the right-hand side of Eq. @ corresponds to Hq. a z a is the Pauli matrix of atom a corresponding 
to the population difference between the excited and ground states, and ojq is the resonance frequency of the dipole 
transition. The second term corresponds to H c g. Here, T — k^d 2 /3ireoh is the spontaneous decay rate of the excited 
state of a single, isolated atom, where fco = 2it/\o = luq/c and d is the dipole matrix element between the ground and 
excited states. <r+ is the atomic raising operator on atom a, and is the lowering operator on atom b. One then 
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sees that the effect of dipole-dipole interactions is an exchange of excitation between pairs of atoms. The strength 
of interaction is modified by a function g(kr ab ) that depends on the distance and orientation between two dipoles. 
It is to be understood that k = fco in the functions / and g. We see that both short-range, near-field (1/r 3 ) and 
long-range, far-field (1/r) dipole interactions are included in our formalism and are treated on equal footing. The 
third term on the right side of Eq. (0) corresponds to L and also is due to atom-atom interactions. It also depends 
on r and has a position dependence described by f(kr ab ). The non-Hcrmitian nature of this term is evident through 
the anticommutator. Physically C describes the processes of both independent and cooperative decay. Finally, we 
have also added phenomenological dephasing through the 7 term, which in particular includes the effects of a finite, 
short-term linewidth of the laser interrogating the clock transition. 

From Eq. J2J, one can derive equations of motion for any atomic operators. For our particular application of Ramsey 
spectroscopy, we find it necessary to solve for the coherence (<x+) and the two-atom correlation (<J^a^). Furthermore, 
to remove the rapid oscillations due to u)q, it is convenient to work in the frame rotating with the interrogating laser 
frequency u>l, where 



dt X a ' 2 x a ' 2 



dt 



= fit + ^±^j (a*a+) T{a+) £ (f(kr ab ) + l g(kr ab )) (*+) Tf(kr ab ){a+al) 



j^a,b j^a.b 

-r E (f(kr aj ) - ig(kr aj )) (a- af of), (5) 

and S — lol — uq. 

In principle, to solve for the atomic system exactly, equations of motion for higher-order correlations are needed, 
although typically some approximation is used to truncate the resulting hierarchy of equations. In general these 
equations can describe light-matter interactions in an optically dense medium, including radiation trapping, level 
shifts, and superradiance 
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III. RAMSEY SPECTROSCOPY IN OPTICAL LATTICE CLOCK 

A. Basic principles 

We now analyze the effects of atom-atom interactions on Ramsey spectroscopy. Starting from the ground state 
i)® N of the system, suppc 
frame by the Hamiltonian 



\g)® N of the system, suppose that one applies a strong probe pulse with the interrogating laser, given in the rotating 



H = E^« Ke* k r " - a-e-*-) , (6) 

a 

where is the Rabi frequency. For simplicity, we have made a plane- wave assumption about the probing laser, taking 
k to be in the positive x direction. We also assume that fewfcoi an d suppress the subscript in future calculations. 
We can do this if the phase error e tSL ^ c over the length of the sample incurred by making this assumption is small. 
Applying this pulse for a time r evolves the system through the unitary operator 



cos Qt e tkx a sm fi T 
-e~ lkXa sin fir cos fir 



(7) 



The state vector immediately following this pulse is given by 

\i/ji)=Y[ (cos^T| ff ) + e ikXa sin flr\e}) , (8) 



where |e) denotes the excited state of the clock transition. 
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In Ramsey spectroscopy, one lets \ipi) evolve for an interrogation time t to p(t). During this time the coherence 
between the ground and excited states acquires some time-dependent phase that depends on 6. After a time t, one then 
applies a second pulse corresponding to the inverse unitary operation W , and then measures the signal corresponding 
to S = ^2 a &a, averaged over the final system pf. S corresponds to the total population inversion. Because of the 
second pulse, S will now depend on 6 and t, allowing one to extract information about the resonance line. 

Formally, we can rewrite S as 

§ = nfcoiuipwuj 

= cos2ftr(^0 -2sin2ftTRe(^y fc:l; "CT+), (9) 

a a 

where the averages denoted above apply to p(t), the system immediately before the second pulse. 
In the case of N non-interacting, independently decaying atoms, 

(o%) = -l + e- rt (l-cos2rh-), (10) 

(e ikx °a+) = isin2fir e -^+ r / 2 +^ 2 ) t , (11) 

which gives a corresponding signal 

S = —N (cos 20r (l - e~ r *) + e~ rt cos 2 2VLt + e -( r +f)*/2 sin 2 20r CQS §t \ _ ^ 

One can see that there is a peak in the signal around 6 = 0. Determination of this peak allows one to find the 
frequency of the transition. One can note two important points about S. The contrast in S with respect to 5 is 
maximized when a "perfect" n/2 pulse is applied, i.e., when Qt = 7r/4. Furthermore, the contribution to S due 
to (cr*) in Eq. © is independent of 8 and thus plays no role in determination of the resonance line. Thus, one is 
motivated to define an effective signal S that consists of the part of S that is actually used to determine the line: 

S= -2sin2f7rRe(^e ifcx »a+>. (13) 

a 

The equation above states that from a theoretical standpoint, determination of the resonance line by measuring the 
population inversion after the second Ramsey pulse is equivalent to measuring the real part of (e lkXa a£) directly 
before the second pulse. 



B. Effect of interactions 



Solving for S exactly in the presence of dipole-dipole interactions appears to be quite a difficult task. Since all 
interactions are proportional to T, our approach is to solve for S as a perturbative expansion in T. In particular, we 
solve for the coherence (<r+) in Eq. to second order in F. This requires solving {(J^<j^) in Eq. (J5J to first order in 
r. The solution for the coherence with appropriate initial conditions is 



where 



B, 



ab 



sin(2£!r)e 



-(i5+ 7 /2)t 



l-y(l + C a ) + ^(l + C7 ) 



+ (/( fcr ^) - i9(kr ab )) (AablH 2 - 2B ab (e~^ +<yt- l)) 



b^a 



A ab = 



|e<*<».-*») cos(2ttr) + I e <*(«.-*») + I (f(kr ab ) + ig{kv a b)) - \f{kv ab ) cos^r) 
t E (/(fcr b ,)-^(fcr b ,))e^^)cos 2 (2f7r), 



j^a,b 



sin 2 (2fir) ^2 (f( kr aj) cos kx aj + g(kr aj ) sin kx aj ) 



j^a.b 



Ca = J2(f( kr *b)-ig(kr a b))e lk( - x ^cos2nT. 

byta 



(14) 



(15) 
(16) 
(17) 
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Eq. I|14fl is correct to every order of 7. Eqs. (|13fl and (|14f> can be evaluated numerically for a given lattice configuration 
and number of atoms. To illustrate the general features of the shift, however, we now make the following simplifications. 
We expand Eq. 114|l to lowest order in 7. We also assume that the Ramsey pulses are nearly perfect 7r/2 pulses, i.e., 
cos2i7r = e<Sl. We then keep terms like eTt but ignore terms like eT 2 t 2 . With these simplifications, 



where 



«>« ^ sin(2fir)e— I 1 - _ + - O ] , (18) 



^ = E (/( fcr ^) ~ ? .9( fcr a 6 )) e 1 *^'-"^ cos(2fir) + (1 + 1 (/(fcr o6 ) + ?5 (fcr ab )) e -**(*.-») 



b^a 



2 X! (/( fcr «j) cosfc;E aj ■+ g(kr a j) smkx aj ) j 



(19) 



The first three terms in the parentheses of Eq. I|18|) are a result of expanding the e _r */ 2 term that appears in the 
result for independent atoms, given in Eq. fT^I . This is just the decay of the signal one would get from independent 
spontaneous emission. The last term in the parentheses is a correction due to atom-atom interactions. Plugging this 
result into Eq. ifH3|) . we find that 



- sin 2 2flr ( (cos St) ( N - ^ + N ^ - E Re^ a j - (sin 5t) E lm<j) a J • (20) 



Because of the antisymmetric sin^t term now appearing in S, one immediately sees that dipole-dipole interactions 
introduce a shift 8 p in the Ramsey fringes, which can be found by solving dS/dS = 0. Suppose that the inequalities 
<5pt<C 1 , rt<C 1 are satisfied. Under these conditions, a simple expression for S p results: 



^ E E (9( kr ab) cos kxab - f(kr ab ) sin kXab) 




1 + 9 E (f( k 

r a j) cos fcxaj + g(kr a j) sin 

J¥a.6 / / 



(21) 



C. Interpretation of shift 

The shift given by Eq. (|21|l yields a simple interpretation. In anticipation of future analysis, we write 8 P as 

(y+T f -)- (22 » 

where 

g(kr ab ) cos fcrc ab - f(kr ab ) sin fcr a &, (23) 

cos20t, (24) 

r a j) cos fcx aj + g(kr a j) sin fciEa?) . (25) 

We will see that U a b is a dimensionless quantity corresponding to the classical interaction energy between two os- 
cillating dipoles, e is a parameter characterizing the error in the Ramsey pulses, and T a is a dimensionless quantity 
characterizing cooperative decay of the system. 

To show the meaning of the U a b term in Eq. I|22|l . consider the interaction between a classical, oscillating dipole 
at r a excited with phase e l ( kx ^-^t) anc j the field incident on it due to a classical, oscillating dipole at v b excited 
with phase e l ( kx b-^t) _ We assume that both dipoles are oriented along the z-axis and that their magnitudes d are 



U ab = 
e = 
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determined from the relation T = k^d 2 /ine^Ti. The classical interaction energy between dipole a and the incident 
field is given by U ab = -(l/2)Re(d Q -E£(r a )). The field at r a due to dipole b is [l^: 

EM . e *..->-g e - (-£! + (W, _ _ , (26) 

where r = |r^ — r a |. Now using the definitions in Eq. J3J), the interaction energy can readily be rewritten as: 

1 TiT 

Uab = --Re(d a • E£(r„)) = — (g{kr ab ) coskx ab - f (kr ab ) sin kx a b) (27) 
= T /7 afc . (28) 

One then sees that this indeed corresponds to the first term of Eq. (1221 . 

Although the U a b term in Eq. I)22ll resembles a classical interaction energy, the terms in parentheses reflect the 
quantum mechanical nature of the system. There is a contribution to the shift that is zeroth order in the interrogation 
time t and proportional to e. One notes that for perfect tt/2 Ramsey pulses, e = 0. Thus, the zeroth order shift is due 
to error in the Ramsey pulse. This can be understood by considering Eq. Q. One sees that the interaction terms 
only influence evolution of the coherence through the term (cr^at). For a perfect tt/2 pulse, this term is initially 
zero, and in this case, interactions cannot affect the measurement at short times. This effect is due to the nature of 
dipole-dipole interactions: these interactions cannot influence the coherence (<r^~) of an atom when it is in an equal 
superposition of the ground and excited states. 

Even if a perfect n/2 pulse is applied, there is an additional contribution to the shift that is first order in t and whose 
strength is given by T a . The intuition behind this is also straightforward. Even if (u^a^) is initially zero, decay of the 
excited state will eventually evolve (a*) away from zero and back towards its equilibrium value of —1. Once (c^a^) is 
nonzero, interactions can influence evolution of (cr+). The rate of decay is characterized by r a . The first term in T a 
is the contribution from independent decay of the atom back to the ground state. The second contribution involves a 
sum over other atoms and represents a correction due to the fact that the decay process may in fact be cooperative 
(e.g., superradiance) . One can easily verify that the contribution from atom j is proportional to Im(d a • E*(r a )): 

Im(d a • E*(r Q ))oc/(fcr aj ) coa(kx a j) + g{kr aj ) sin(kx aj ). (29) 

This reflects the well-known result that the atomic inversion (<t„) is driven by the dipole component in quadrature 
with the incident field. 



IV. GENERALIZATION OF RESULTS 



A. Imperfect filling of lattice sites 

Experimentally, knowing the exact number of atoms in the lattice and achieving a filling factor of one atom per 
lattice site are difficult tasks. Most likely, one can experimentally determine the density p(r) of atoms in the lattice, 
such that the probability of occupation at any particular site a is P{v a ) = p{Y a )V 1 where V is the volume of a unit 
cell. It is straightforward to modify Eq. I|21|) to the case of imperfect filling. For simplicity, we only consider the shift 
that is zeroth order in t. This shift can be written 



Y = ^^2^2^cos2nT(g(kr ab )coskx ab - f{kr ab )sinkx ab ) (30) 

a b ^a 

2 1 

= — -cos2nr g(kr ab )coskx ab (31) 



N t-r 1 2 

pairs 

= I C os2r!r^C/(R)^^, (32) 

R^O 

where J7(R) = g(kR) cos kR x , {R} denotes the set of direct lattice vectors, and N(R) is the number of pairs of atoms 
separated by R. In the derivation above we have utilized the fact that sin kx ab — — sin kx ba to cancel the sum of 
f(kr ab ) sin kx ab . In a realistic scenario, one neither knows N(R) nor N exactly. In this case, one must solve instead 
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for the ensemble average (S p ) and the variance AS p . For large N, one can safely pull the factor of N out of the 
ensemble average: 



,AT(R) . <JV(R)) 
1 JV '"^ (AT) 



With this simplification 

r 



-^cos2r!r^[/(R)<7V(R)) 

^ ' R#0 



1- cos2r!T ^ U(R) J dvp(v)p{v + R)V, 



2(N) 



(33) 

(34) 
(35) 



(A5 P ) 5 



R^O 

2 



^cos20rj ^ f/(R)[/(R')(iV(R)iV(R')) 



^_cos20r53f/(R)(iV(R))j 



(36) 



1 



2 (TV) 



cos 2Qt 



53 C/(R)C7(R') /drp(r)p(r + R)y 2 ((l-p(r)y)p(r + R') 

R,R'^0 "* 

+ (1 - p(r + R)V)p(r + R + R) + (1 - p(r)V)p(r - R') 
+ (1 - p(r + R)V)p(r + R - R ) 



+2 53 U{R) 2 J dvp{v)p{r + R)v((l - p(r)p(r + R)F 2 ) - (1 - p(r)V>(r + R)V 



R^O 

-p(r)(l-p(r + R)VOV 



(37) 



Each term in Eq. 1|37|) has a clear meaning. To calculate the average shift in Eq. I|35|l . one must evaluate (iV(R)). To 
do this, one must perform a sum over of the probability that the sites and + R are both occupied. To find the 
variance, one must calculate quantities like (AT(R)AT(R')) , and thus the probability that the sites r,, + R, r,-, and 
Tj + R' are all occupied. When these four points are distinct, the probability is simply a product of the probabilities 
of each point being occupied. This is untrue when one or more of the points overlap. The terms in Eq. I|37|l represent 
corrections due to these overlaps. The product p(r)p(r + R)V 2 ((l — /j(r)V)/o(r + R'), for example, is due to the overlap 
of r, and r, . 



B. Effects of multilevel atomic structure 



Our results derived thus far are for the case of two-level atoms. This is the relevant case of study for the (J = 
0) — ( J = 0) forbidden transition proposed for optical lattice clocks, where the simple level structure makes it easier to 
cancel the relative Stark shift in the clock transition. Nonetheless, our results can be generalized to more complicated 
level structure, such as the case of an atom with a single ground and multiple excited states. A simple argument 
shows that Eq. I|21|) remains correct to the lowest nontrivial order in Ft. If multiple excited states are present in 
addition to the one that is initially excited, the equation of motion for (erj clock) i n Eq. 1|18|) will contain additional 
terms like (cr*0^~ o thcr): where the subscript "clock" refers to the clock transition and "other" refers to other excited 
state levels. Initially, (another) = and thus 

(pother) ^ (3g) 
{V^b dock) 

Consequently, at short times, evolution of (cr^~ c iock) will be dominated by the clock transition. Thus, if imperfections 
in the Ramsey pulse constitute the major source of shift, the multiple excited states will contribute an additional 
source of shift that is first order in Ft. If decay of the clock excited state constitutes the major source, the multiple 
excited states will contribute a shift that is of order T 2 t 2 . 
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V. ANALYSIS OF RESULTS 



Eq. H21|) or Eqs. (|35[l and (|37|l can be evaluated numerically for a given lattice configuration and number of atoms. 
To extract the key features of the shift, we note that the zeroth-order shift in t in Eq. I|21|) essentially consists of adding 
together the classical dipole interaction energies f/ a fcOc — Re(d Q -EJ(r a )). For a generic configuration of atoms, the 
amplitudes of the dipole fields incident on a given dipole a tend to interfere. For certain configurations, it is possible 
that the field amplitudes will add constructively along some direction of propagation k. Near these configurations one 
will expect large shifts to result. The condition for constructive interference between radiated dipole fields is similar 
to that of Bragg scattering in a crystal, and is readily found to occur when 

|G|=fc , (39) 
where G = (G x — k , G y , G z ) and G is a reciprocal lattice vector. This condition can be rewritten as 

|G| 2 = 2k G x . (40) 

Numerical results indicate that peaks in the line shift do indeed occur when condition H4U|) is nearly satisfied. 

One can easily derive an approximate scaling law for the line shift in these resonant configurations. We define 
a dimensionless parameter (3 related to the density of atoms by n — l/(/3A) 3 . /3 characterizes the spacing between 
neighbors in the lattice. In a resonant configuration, the electric fields add constructively, and the total electric field 
experienced by an atom is approximately 

f , n f L nr L 2 

where L is the linear size of the system. For N total atoms, L ~ (3XN 1 / 3 . Then 

S N 2/3 

^ ~ (42) 

Experimentally, one has freedom to choose the orientations of the trapping laser beams that form the lattice. The 
control parameters can be parameterized by a set of variables {a}, which will also determine the reciprocal lattice 
vectors G({a}). One can then find solutions {cto} of Eq. I|4(J|I corresponding to configurations with large line shifts. 
In the parameter space between two sets of solutions {ao}; one can numerically find configurations where the shift is 
significantly reduced. 

In the case of imperfect filling of lattice sites, it will be important to account for not only the mean shift but the 
variance as well. For large numbers of atoms, Eq. 1)37(1 cannot be evaluated exactly without extensive computational 
resources. With a small filling factor P — pV^ll, however, we can estimate that the major contribution to the 
variance results from the p 2 terms, while the p 3 terms remain negligible. In this diffuse limit, 



(AS p ) 2 r 



cos 2^ x2 ^(R-) 2 / drp(r)p(r + R)V. (43) 



r 2 \2{N) 

In this case, one readily finds that the variance A5 P scales like (P/N) 1 / 3 . For P < 1/2, the variance increases with 
P due to the increasing uncertainty of whether a pair of sites will both be occupied, but decreases with N due to the 
decreasing fractional uncertainty in the total number of pairs of atoms separated by a vector R. 



VI. NUMERICAL EXAMPLE 



As an illustration of our results, we consider 87 Sr atoms trapped in a lattice formed by six interfering beams, as 
shown in Fig.^ For 87 Sr, the "magic" wavelength of the trapping lasers is roughly Xl — 1.07Ao Q, and one can vary 
the angle 6 between the propagation vectors of the trapping beams. The resulting lattice is tetragonal, with lattice 
constants of a x = ir/kLsmO, a y = n/kLCOs8, and a z — Tr/fc^ along x, y, and z, respectively. The lattice constants 
are plotted in Fig. The corresponding basis of the reciprocal lattice has lengths G x — 2fcisin6', G y = 2A;_lcos#, 
and G z — 2k^. We have ignored the effect of atomic back-action [l3L IT^j on the trapping fields, whereby scattering 
of light by the atoms introduces phases that might modify the lattice constants. Such effects are expected to be 
stronger in red-detuned lattices, where atoms lie in the antinodes of the potential, and with increasing atomic density. 
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Taking into account this back-action does not modify our results, except that now the lattice constants must be solved 
self-consistently OH- 

Using Eq. (|4C)|1 we can find values of 9 where constructive interference causes the shifts to be large. We focus on two 
specific solutions, 9o/ir = 0.116 and 0o/tt = 0.180. For our system we consider (N) atoms in a spherical distribution 
with uniform density p(r) for r < ro, and zero density for r > r^. The relationship between the density and filling 
fraction P is given by P — p(r)V, where V is the volume of a unit cell. The critical value ro is determined by the 
equation 

ro= [-^p-) ■ (44) 

We first consider a perfectly filled lattice consisting of N = 10 6 atoms. For simplicity we calculate the line shift to 
zeroth order in the interrogation time t. In Fig. ©, we plot the quantity 25 P /T cos(2f2"r) as a function of 9. Peaks in 
the shift are clearly visible at the points 9q that were calculated analytically. It should be noted that the line shift 
can be very large in one of these "bad" configurations. Even in the limit of short interrogation times, one can see 
that shifts of order (5 p ~10r are possible. This is perhaps a surprising result, and occurs because the spatial ordering 
of the atoms allows the interactions to behave constructively at these points. For longer interrogation times, one 
expects this line shift to become even larger, since the constructive interference in these configurations also leads to 
superradiant decay and thus a large contribution to the shift that is first order in t. One also sees that away from 
these "bad" points, the shift is strongly suppressed and even becomes zero for one particular value of 9. 

We next consider a partially filled lattice consisting of (N) = 10 5 atoms and a filling factor of P = 0.05. Fig. 0] 
gives the quantity 2(6 P )/T cos(2J7r) as a function of 9. The mean shift (6 P ) still exhibits peaks at the points 9 , and 
vanishes near 9/w = 0.125. At this "good" point, one can use Eq. (|43|l to estimate the variance in the expected shift. 
Within this diffuse approximation, we find that the variance 

^«icos(2ftT)x3.1xlO~ 3 . (45) 

Experimentally, there will be additional sources of error that result from not knowing p(r) perfectly, errors in the 
configuration of the trapping lasers, and from the effects of atomic back-action on the lattice constants. Nonetheless, 
it appears that the error due to dipole-dipole interactions can be made quite small by appropriately designing the 
lattice. 



VII. CONCLUSION 



We have derived an expression for the line shift measured in Ramsey spectroscopy due to dipole-dipole interactions. 
We find that the lattice geometry strongly affects the magnitude of the shift, and is peaked in lattice configurations 
where the interactions between atoms add constructively. Because of the spatial ordering in the lattice, the shift can 
be quite large in these resonant configurations. By tuning the lattice between two of these configurations, one can 
reduce the dipole-induced line shift to nearly zero. 

While the resonant configurations might be bad for clock applications, it might be worthwhile to study these 
configurations further. The dipole-dipole couplings in an optical lattice offer the possibility of strong, constructive 
interactions that can be dynamically tuned by changing the lattice geometry. This might be useful for applications 
such as quantum information processing and might have interesting consequences for studying phenomena such as 
superradiance and for probing the superfiuid-Mott insulator transition. 
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APPENDIX A: DERIVATION OF MASTER EQUATION 



In this appendix we derive Eq. © starting from the full atom-field Hamiltonian. For a more detailed derivation 
and discussion, one can also see [l5l Il6| . The Hamiltonian for the atom-field system is 

H = H + V, (Al) 
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where 



Ho = Hi 



internal ~l~ Afield 



k. c : 



and 



v = -E d r%) 



(A2) 
(A3) 



(A4) 



(A5) 



cr| refers to the internal state of atom j, luq is the atomic transition frequency, and Wk.e is the frequency of the mode 
of the electromagnetic field with wave vector k and polarization e. Tj is the position of atom j, and e a t om is the 
polarization direction of the dipoles. 

We now consider the evolution of the atoms+field density matrix p a f in the interaction picture. The equation of 
motion is given by 



dt 



1 

it) 



V(t), Paf 



V(t) — e iH ° t / h y e ~ iH ° t / h 



(A6) 
(A7) 



We can integrate Eq. <| A6|> once and substitute the result back into itself. We then trace out the field degrees of 
freedom to obtain an equation of motion for the atomic density matrix p alone: 



dp 
di 



1 f 
— Tr/ / dr 

h Jq 



V(t), V(t-r),p af (t-T) 



(A8) 



To make the above equation more useful, we employ the Born-Markov approximation, replacing p a f{t — r) with 
p(t) (g> |0)(0|. Physically, this amounts to assuming that correlations in the field are negligible, that the field can 
always be approximated by a vacuum state, and that the correlation time of the atom-field system is much shorter 
than any atomic timescales. These assumptions safely allow us to extend the time integral to infinity, so that 



^ = 2^ r / / dl ~ 

h Jo 



Of 



V(t), V(t-T) )P (t)®\0)(0\ 



Writing out V(t) and V(t'), where t' = t — r, give 

V(t) = -^k,e (eatom • e) (ir+e^* + u~ e'^ 

a, k,e 

v(t') = J2 • e ") (°t^ ot ' + °b^ ot ' 

b, k',e' 



ak, e e i(k ' r "" a ' k '« t) + 4 £ e- i(k - r »~ Wk ' et; 



(A9) 



(A10) 



(All) 

When we substitute the expansions above into Eq. i|A9[l and perform the trace, the only nonzero terms will be those 

associated with ak, e 4' e / 10) (0| , 4 E |0)(0|ak', e '> 4' e '|0)(0| a k,e, and |0)(0|ak', e '4 c wnere k = k' and e = e'. Making 
these simplifications and replacing the sum on k by an integral gives 



dp i r°° 



(2tt) 3 



k 2 dk dnJ2d 2 £lA^tora-£) 



e -i( a)k(e - Wo )r+ik.r oi , a .+ cr - p 



i „-»Kt+"o)T+ik'r a i.-.+ . _ „!(w k ,,+wo)T- 
t e a a a h p e 

_ e i( Uk ..-<-o)T-ik.r a » (7 .- p(T + + h( . 



£k-r a 



(A12) 
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The angular integral is tedious but straightforward and gives 

/\ V^ ik-r ab t~ ^2 . ( . 2a smkr ab , 2a \f coskr ab sinfcr Qb \\ 

J dn ? 6 (£atom ' e) = 47r l sm + (1 - 3 cos 0) It^f " Jk^r ) ) ■ (A13) 

The time integral can be evaluated using the formula 

f°° 1 

/ cdre- lc(fe=Ffco)T =Tt5(k Tk)±iV- r , (A14) 

Jo k oTk 

where V denotes the principal value. The delta function above eventually yields the non-Hermitian component of the 
evolution, while the principal value yields the coherent atom-atom interactions. Performing all integrals, making the 
replacement T — k^d 2 /3ireo?ic 3 , and changing back to the Schrodinger picture yields Eq. J5J). 



APPENDIX B: SECOND QUANTIZATION OF ATOMS IN LATTICE 

The results above were derived treating the atoms in the lattice in first quantization. Starting instead from second 
quantization, one can sec that the results derived in first quantization are appropriate in the limit of tight confinement 
of the atoms, when the overlap between atoms at different sites can be ignored. 

In second quantization, the atomic wave function ^>(r) can be expanded in terms of Wannier functions: 

^( r ) = ^2 b vni<f>vn{r- r i)- ( B1 ) 

vni 

Here v = e,g denotes the internal state of the atom, n is the band index, and i labels the lattice sites r,. To lowest 
order, we assume that only the harmonic oscillator ground-state wave function is relevant. The atomic Hamiltonian 
is then given in second quantization by 

vij 

In the limit of tight confinement, the overlap integrals for i^j can be ignored, leading back to the atomic Hamiltonian 
used in first quantization. 

In second quantization, the electric dipole Hamiltonian is given by 



H ' = ~ E ( b ie b 09 J dr Mr - ri)(d • E(r))^(r - r,) + b\ g b je J dv^v - Fi )(d • E(r))0 e (r - r 3 -)) , 



(B3) 



where 



E(r) = ^k,ea k , e e* r e + h.c. (B4) 

Again, in the limit of tight confinement, the overlap integrals for i^j can be ignored, and furthermore the term 
E(r) can be replaced with E(r»). In this limit this Hamiltonian is equivalent to the electric dipole Hamiltonian used 
in first quantization. 
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FIG. 1: The lattice studied in our numerical example is formed by the interference of six laser beams. Thick arrows denote the 
directions of propagation of the beams, and thin arrows denote the direction of polarization. Four beams are oriented along 
the x-y plane, each making an angle ±# with the t/-axis and polarized along z. Two additional beams run parallel to z and are 
polarized along x. 




0.1 0.2 0.3 0.4 0.5 

e/7t 

FIG. 2: The lattice constants for the six-beam lattice are plotted in units of Ao as functions of 9. We assume that the ratio of 
the resonant wavevector to that of the trapping lasers is ko/kL = 1.07, consistent with the magic wavelength of 87 Sr. The solid 
line represents the lattice constant along x, the dashed line along y, and the constant dotted line along z. 



14 




FIG. 3: The calculated shifts for the six-beam lattice as a function of 9. The system consists of N = 10 atoms with a filling 
factor of 1. 




FIG. 4: The calculated mean shifts for the six-beam lattice as a function of 8. The system consists of (N) = 10 5 atoms with a 
filling factor of 0.05. 



